PRODUCT FORMULA FOR JACOBI POLYNOMIALS, 
SPHERICAL HARMONICS AND GENERALIZED BESSEL 
FUNCTION OF DIHEDRAL TYPE 



NIZAR DEMnA 



Abstract. We work out the expression of the generaUzed Bessel function 
of i?2-type derived in This is done using Dijskma and Koornwinder's 

product formula for Jacobi polynomials and the obtained expression is given by 
multiple integrals involving only a normalized modified Bessel function and two 
symmetric Beta distributions. We think of that expression as the major step 
toward the explicit expression of the Dunkl's intertwining operator Vk in the 
B2- invariant setting. Finally, we give in the same setting an explicit formula 
for the action of Vk on a product of lyp^jK > and the ordinary spherical 
harmonic KimCy) := Ij/I*"" cos(4m(9), y = |j/|e*^. The obtained formula extends 
to all dihedral systems and it improves the one derived in |16| . 
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1. Motivation: Dunkl's intertwining operator 

The reader is referred to [7j and references therein for an extensive background 
on the theory of Dunkl operators. The most challenging problem in this theory 
is an explicit expression for the action of the so-called Dunkl intertwining opera- 
tor denoted in literature by Vk- The latter is a linear isomorphism of the space 
polynomials in several variables which is degree-preserving, conservative and which 
intertwines the commutative algebra of Dunkl operators and the one of partial 
derivatives. Few, yet relevant, works aiming to solve the above-mentioned problem 
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were achieved and complicated formulae for the action of Vk were supphed ([7], [S], 
[17j ) except in the orthogonal case corresponding to products of copies of the group 
Z/2Z. In that case, Vk is represented by means of a multivariate Beta distribution 
which reduces to a symmetric one (invariant under sign changes maps) in the W- 
invariant setting {W is the reflections group), while for other nonorthogonal root 
systems multiple integrals with complicated integrands were obtained. Neverthe- 
less, a relatively easy and elegant formula, compared to others, for the action of Vk 
on (x,-), where a; G is fixed and (•, •) is the Euclidean scalar product, was given 
in [in] for the root system of type i?2- This root system is the most elementary 
example of non orthogonal root system and matches with the dihedral system of 
order 4, denoted by /2(4) and corresponding to the symmetry group of the square. 
A noticeable coincidence is that while the proof supplied in [16^ relies heavily on 
the product formula for Jacobi polynomials due to Dijskma and Koornwinder ([5]), 
the expression derived in [4J for the generalized Bessel function of dihedral type, 
say , of dihedral type also involves a product of Jacobi polynomials. Thus it 
is quite natural and interesting to use the product formula for Jacobi polynomials 
in order to seek a more elegant expression for in the -B2-type setting, aiming 
to investigate the connection between both works [4J and 6^ and to have a better 
insight into the action of Vk on _B2-type invariant functions. In order to motivate 
the reader, we show below how the easy product formulae for cosine and sine func- 
tions allow, in some particular cases, to get an easy expression for of -B2-type. 
The general setting, corresponding to general Jacobi polynomials, is more subtle 
and one is far from getting (even from hoping) an easy expression for as in 
those particular cases. Nevertheless, the formula we supply involves only a modified 
Bessel function and two symmetric Beta distributions and simplifies considerably 
for some particular values of the multiplicity function, yielding as a by-product to 
the positivity of . 

1.1. Particular cases. Recall that even dihedral groups D2{2p),p > 2 are the 
symmetry group of a 2p-gone and that it contains two classes of reflections, with 
mirrors joining opposite vertices and mirrors joining the midpoints of opposite sides 
([?]). The following formula was derived in |3|: let a; = pe'*^, y = re'^ where p,r > 
and < 0, 61 < 7r/(2p), then 



where k = (ko, ki) is a positive- valued multiplicity function, U — ki — 1/2, i g {1, 2} 
are the index values, 7 = p(fco + ^i) and IiJ,Pj' " are the modified Bessel function 
and the j-th (orthonormal) Jacobi polynomial respectively. A similar formula holds 
for odd dihedral systems D2 (n) , n > 1 where one has to substitute in (|l.ip ki = 
0,p = n, fco = ^ > 0. The constant Cp^k is such that 13^(0,?;) = 1 for all y and 
is easily computed as follows: recall the series expansion of the modified Bessel 
function: 



(1.1) 

Df{x,y)^DY{p,cj^,r,e) 




J2 hjp+, {pr)p'l (cos(2p0))py (cos(2p0)) 



3>0 



(1.2) 




2 



Thus, for p — 0, the only non zero term in the double infinite series (after substi- 
tution of <\l-'2\i ) corresponds to j = m = 0, therefore 

c;l = ^;^p[,-'°(cos(2p0))p[,-'°(cos(2p0)). 

Since Pq is constant and has unit-norm in 1, 1], (1 — + uf"du)^ then 
\ _ oko+kiQQ^^ + 1/2, /co + 1/2) 



p[,^''«(cos(2p0))p[,-'°(cos(2p^)) 



so that 

ri3l c - 2^"+^^ ^^P^^^ + ^"^ + + l/2)r(fco + 1/2) 

^ ' ' "'^ r(A:o + fci + l) 

The particular cases we mentioned above correspond to fco = = and fco = 
ki = 1 and are referred to as the geometric cases. In the former, reduces 
to the orthonormal Tchebichef polynomial of the first kind defined by (|llj) 

T,(cos0) = y|cos(j0), j ^ O,ro(cos0) - ^ 
so that I?^ simplifies to 

(1.4) Z?o^(p,0,r,0) = ^ 



/o(/5r) + 2 ^ hjpipr) cos{2jp<j)) cos(2jp6') 

= A)(P?')+ hjpipr) cos{2jp4>) cos{2jpe) 

jez\{o} 

= X! hop^P"^) cos{2 jpc/)) cos{2 jpO). 

since Cp :— Cp^ — tt and where the second equality follows from /_„i(z) — Im{z) for 
m e Z (which is a consequence of J-m{z) = (— l)'"Jm(z) and Iu{z) — {—lY Ji,{iz), 
see [TS], p. 15). Using the product formula 2 cos a cos & = cos(a + fe) + cos(a — fe), one 
is led to compute (both series converge as we shall see) 

E I^jpipr) cos(2jp(0 ±0))^IJ2 [e»2,p(^±«) ^ g-.2,K0±«)] ^ ^ ^^^.^(^^)^ 

Now, let ^ = e^^^P and j G Z, then it is an easy exercise to prove that 



S — 1 



otherwise, 



for every p >2. It follows that 

E^2,,(pr)e--(^±«) = lf;^/,.(He^(^±^)^-r 

jGZ ■'^ s=l j£Z 

Now using the generating function for Bessel functions (fT^ p.l4) 

(1.6) e(*-V*)^/2 = ^J,(z)^^ teC\{0}, 

3 



together with Ij{z) — {—iyjj{iz), one gets 

-, 2p 2p 

jez ^ s=i jez ^ s=i 



pr cos[(0±^) + s7r/p] 



As a result 

gpr cos[(0+e)+s7r/p] ^ ^ 



^prcos[{<j>-e)+STv/p] 



_s=l s=l 

which is in agreement with the definition 

and with the fact that the dihedral group contains 2p rotations (a; i~> xS,'^ , 1 < 
s < 2p) and 2p reflections (a; i-^ x^'', , 1 < s < 2p). Note that for p = 2, since 
s7r/2 e {7r/2, TT, 37r/2, 27r}, then easy computations transform D^^ to 

a.7) 

- [cosh(pr cos(0 + 6)) + cosh(pr cos(0 — 6)) + cosh(pr sm{(f) + 9)) + cosh(pr s'm{(j) — 9))] . 

Similarly, in the latter (fco — ki — 1), p^^'^^^ reduces to the (orthonormal) 
Tchebichef polynomial of the second kind Uj defined by ([2]) 

t/.(»s(.)) = Ji:i^M,.>o. 

V TT smb* 

Using the product formula 2 sin a sin & — cos(a — &) — cos(a + 6), one finds with 
u;(r, 9) r^P sin(2p6l) that 

2p 

y 

The last expression agrees with 



}_ \ ^ I prcos[{rf>-9)+STT/p)] _ pr cos[{rf>+6} + STr /p)] 



u;{r,9)u;{p,4') ^ 

which follows from the shift principle ([?]) or the reflection principle and the Doob's 
transform ([Sj, [l2]). 



1.2. Results. For general multiplicity values, the situation is far from being easy 
since for instance, of type B is given by a multivariate series qFi defined 
via Jack polynomials ([5]) whose explicit expressions are not easy to write down 
(though note the recent progress in this direction) except in the few geometric cases. 
Nevertheless, we succeeded to express of i32-type as a multiple integral involv- 
ing only a normalized modified Bessel function and two symmetric Beta probability 
measures depending on ki, kg, that is we got rid of the dependence of on Jacobi 
polynomials. More explicitly, let 



, 1 + u COS 2(? COS 20 + V sin 20 sin 20 , 
C20,2e(w, v) := W ^ , [u, v) e [-1, IJ , 

and, for t e [0,1],Q < z < 1/2, let 

1 /"l g7/2-l ^ 

Kj{t,z) := —= / — =^=«^/2(iVl - 2z(j)d(7, 7 > 0, 
i,y being the normalized modified Bessel functiorQ defined by 

(1-8) ^-(^)^=E w ^ , (o) ■ 

"'-^ 1 m + 1/ + 1 m! \2/ 

m— 

Denote also /i" the symmetric Beta distribution 

(1.9) ^i^du) (x{l-u'^y-'^^'^l[^i^i]{u)du, v>-\l2. 

Then, our first result may be stated as follows: 

Proposition 1.1. 

i^r(x,,)=r(l±i) I / ' + (prc2,,2.(^, V)) + sin^(77r/2) 

(1.10) 



7^7 I P*^*. ^^^^^ j + —:^d,K^ I prt, — ^^^-^ 



dt.fi'^{du)^i'°{dv), 



The extension to j = is performed via weak limit and one recovers (|1.7p . 

Remark. p.lOp is relatively easy and explicit regarding the existing formulae (\8J, 
|17j j. Besides, we think of it as a major step toward an explicit expression for the 
action ofVk on B2-invariant functions since Bessel functions enjoy a huge number 
of nice properties. However we do not know how to come to Vk due mainly to the 
last term in the integrand. 

Corollary 1.2. For even values 0/7 7^ 0, 

D^{x,y)=T(^^^ J y Z(^-i)/2 {prc2^Mu,v))fi''{du)f,'»{dv) 

which is obviously positive. 

Our second result supplies for all dihedral systems, independently from the first 
one, an explicit formula for the action of Vk on a product of > and 

the ordinary spherical harmonic Ym{y) = |?/|™ cos(m6'), ?/ = |2/|e*^. For sake of 
simplicity and coherency with the first result, we only write down the formula for 
the D2(4)-type root system: 

Proposition 1.3. For positive integers K,m and y = \y\e''^ , one has 

Vk [I • rn™] {y) ^ \y?^^'- E r V-vrn'^r-t" t 1^ ^^+- 

(k - 2:;)!r(4m + 2j + K + 7 + 1) ■'^ 

where bm,j+m is the m-th coefficient in the expansion 0/ P^'^''° (cos(4^?)) as a finite 
linear combination of cos{Am9) . 



^This function is slightly different from the spherical modified Bessel function which is valued 



1 at 0. 



Remark. Proposition (II. 3p shows the action of on polynomials that are invari- 
ant under the group i?2- Such results have been studied and can be derived after 
some manipulations from 16], but the formula obtained here is more explicit. We 
didn't attempt to compute the special cases in which the explicit formula ofbj „i 
known. This is of interest only if the formula can be further simplified or the sum 
in the proposition can be summed up. Note also that our formula resembles the one 
displayed in Theorem. 2 p. 13 in [10 . 

2. i32-TYPE ROOT SYSTEM: PROOF OF PROPOSITION (|1.1[) 

As indicated above, this section consists of six subsections, the title of each 
one indicates its content, aiming to orient the reader through the overwhelming 
yet tricky computations. Before coming into them, we want to inform the reader 
that usual operations performed on integrals (exchange of integration's order or 
limiting operations under integral signs) are easily justified due to the compactness 
of integration domains and to the smoothness of the involved functions. 

2.1. Product formula: a first transformation. The generalized Bessel function 
displayed in (jl.ip reads in the i32-type setting 



D]^{x,y) = C2,k(-y ^j+7(pOp;yi''(cos(40))p;.-|''(cos(40)) 



where 



937/2 + y 
= 1^ ) r(fci + i/2)r(fco + 1/2) 



by Gauss duplication's formula. Now, recall the product formula for Jacobi poly- 
nomials ([iJJ [m p.4245) 



c^^ppf^{cos2c^)pf''{cos2e) = {2] + cy + p+l) j j C^+^+\zeiu,v))^J,"{du)^^^{dv) 

where ^{a),^i(3) > -1/2, (j), 9 e [-7r/2, 7r/2], 

_ ,„+g+i r(a + l)r(/3+l) 
''"^^"^ r(a + /5 + l) ' 

(2.1) z^^g{u, v) — u COS 9 COS (j) + V sin sin (p 

and /i" is the symmetric Beta probability measure whose density is given by 

W = 0Fr(a+V2) ^^ ~ «^)"-i/2l[-i.i](u)d^, a > -1/2. 

Speciahzing the product formula to a — li — ki — 1 /2 > —1/2, f3 — Iq — ko — 1/2 > 
— 1/2 gives 



J=0[4] 

r(7)(f)^ E Ij+i{pr){j+l) I I CY+''-{z2^,2eM)^i'-{du)^i'^^{dv) 



^We write the formula using orthonormal Jacobi polynomials and symmetric Beta probability 
measures for later purposes, yet the Gegenbauer polynomial is not normalized. 
"^There is an erratum in the constant term in front of the integral. 



by Gauss's duplication formula. Now, a useful by-product of the product formula 
is the following (Theorem 2.2. in [B] specialized to a = /? = A — 1/2, A > 0) 



C^iz) = J Ci} ^'-^f^^j ,'-"'{dw) 
and yields with A = 7/2 

Dr(x,t/)^r(7)(^^y /,+,(H(j+7) / / Jc]iz)f/Hdu)^.'»idv)^,'^■'''^/'idw) 



i=o[4] 

where we set 



Z = ^20,2e(w, w, w) := Y 2 ^■ 

Finally, using (jl.Sp . one gets: 
(2.2) 

2.2. An auxiliary formula. Since 

(2.3) |c7(Z)|<|C7(l)| = i^, 
and since 

E^^.+.(H <oo, 

then Fubini's Theorem applies and one is led to compute 

^(j+7)/,+,(pr)C;(Z)e^ 

for fixed (u, w, w) e] — 1, Note that for s = 2, 4, one has 

(2.4) ^(±l)^-(j +7)/,-+,(pr)C;(Z) = ^^-f^"'^ 

by formula 5.13.3.3. p. 712 in [2], yet we did not find any similar result to 

Y,{±tYij+-/)I,+,{pr)C]iZ). 

However, note that for strictly positive integer values of 7, (|2.4p may be written as 
J2{±inj + 7)/,+.(pr)C;(Z) = ^^^Dl [e^P^^] . 

In fact, we can derive a more general similar result for 

Y,{3+l)I,+,{pr)C]{Z)e^^'^ 

7 



for any real number q provided that 7 > 1 is an integer (then extend it to strictly 
positive values). In fact, from p. 32 and p. 44 in [T?], one gets 

2^/,(^)e*^^ = 2^/,(^)[cos(X) + ^sin(jC)] 

j>0 j>0 



= exp(zcosC) + -/o(C) + *sinC / exp(— fcosC)-/o(0'^^ 

JO 

for all real numbers z, C and from 11.1.2 (18) p. 235 in [llj, the following holds 

(j + 7)C7(Z) - —^^Dl{T,+,{Z)) = ^_l^i5|[cos((j + 7) arccosZ)], 
where Dz stands for the derivative operator. Moreover by the bound (|2.3p . one has 

E(J' + l)h+,ipr)Cl iZ) e^'^ = Ij+,ipr)Dl[cos{{j + 7) arccosZ)]e^^'' 

= 27-ipc^) -^z ^jipr) cos(j arccosZ)e'J«, 
and since cos(j arccos Z) is a (Tchbeycheff) polynomial of degree j, then 

+ l)Ij+'r(.pr)C] (Z) e^^" = J.'Z. Dl ^ I,{pr) cos(j arccos Z)e'^^' 
j>o j>o 

where we set a := cosZ. 

2.3. Odd and even values of 7. Easy computations using the last above equality 
yield 

J2{±iYU + l)I,+,{pr)C]{Z) = ^^^Dl cosh(prVl - Z^)] ±iZ cosh{W 1 - Z^)h{t)dt 
for s = 1, 3. If 7 is even, then i'^ — so that 

(_l)7/2 



(2.5) i+)J2i3+l)I,+,{pr)C] (Z) i±^y^ = ^' Dl[cosHprVT^)], 

j>0 



while if 7 is odd then {—i)'' = —i(~iy' ^ = ~i{iy' ^ ~ '^i^T' therefore 
(2.6) 

(_l)(7-l)/2 rPr 
J ^l)^3+l\Pn^j \^)\^'^)- = n7r^ . 

J>0 



(+)E(j' + ^)^^+^(Hc^7(^) = ^ 27-ir(7) ' ^"^ti(t yT^)/o (Od^j ■ 



Above we used the symbol (+) to indicate that we sum both series. For s = 2,4, 
one recovers (|2.4p and both series displayed there contribute to 



{+)Y,{±lY{3+^)I,+,{pr)C]{Z) = J^'l' cosh(prZ). 



Let in denote the normalized modified Bessel function already defined in (jl.Sp . 
it has the Poisson integral representation for v > —1/2 ([TT] p. 81): 

z,(z) ^ ^ f e^-(l - w'y-^/'dw = — ^ /' cosh(z«;)/i-(d«;). 

Therefore, one gets with = (7 — l)/2, 
(2.7) 

(pr) 



7 



|^^(+)^(±1)^"0- + 7)/,+^(pr)C7(Z)A.(^-^)/^(d«;) = 2r ^(,_l)/2[prc20^2.(^^, «)], 



where we set 



„ /I + 22,^,29 (W,?^) , X 

Z = \j ^ w := C2^.2b[u, v)w. 



With regard to (|2.5p . one needs then to integrate 

(2.8) Dl[cosh{prVT^)]^z=c.^Mu,v)..f^^^''^^Hdw), 

for even 7 > 2. 

2.4. A positive-definite function. One tricky way to compute (|2.8p is to use 

the positive-definiteness of Z cosh(t-\/l — Z"^) for all real t. In fact, the following 
Bochner representation holds and is easily derived from 6.645.3 in [13^ by an analytic 
continuation 



cosh(iyT^Z^) = cos(tZ) [' e'tZq hi Wl q ^^^ _ j 
where vt is the symmetric measure 



2 7-1 



^ ^ Vi - 

Thus, one has for even 7 

cosh(prVl - = (-1)7/2(^^)7 j q^e''^^vpr{dq). 

Using Fubini's Theorem, it follows that 
»i 



(_l)7/2 j ^ [,osh(prx/T^)]|2^,,^ ,,(„,.)„^(-^-i)/2(dw;) = 

(pr)'^ J q^ (^j e*'"^=^*-^«("'")"'V^''"^^^^('^^«)^ i^pr(dg) 

and the integral between brackets is nothing but the spherical Bessel function 
J(7-i)/2 defined by j^^z) :— in{iz) so that: 



7 + 1 



2 

9^i(7- 1)/2 [C20,20 (m, W)p?-g] i^pr (rf?) ■ 



With regard to (^3]) and pTT)) . 

2^r(7) 



(pr)T 
(2.9) 



7+1 



Now, using the expansions of j^_i/2 and /i, one has 

»i 



y 



-1 Vl-?" j,m>0 



<f j{i-i)/2{c24>,2e{u,v)prq)vpr{dq). 
(-ly 1 



r(j + (7 + l)/2).?! (m + l)!m! 



E 

J , m > 



r(j + (7 - l)/2 + m + 2)j\ (m + 1)! 



1 / C24,^20{u,v)pr\'^' ( pr 



E' 

i>0 



(-1)^ / C24,,2eiu,v)pr 



2j 



E 

m>l 



r(j + (7 - l)/2 + m + l)m! 



(f) 
(?) 



2m+2 



(-!)■' / C24,ae{u,v)pr 



2j 



pr 



(7-l)/2+J 



-'j + (7-l)/ 



2(P?') 



pr 



(7-l)/2 



r((7-i)/2 + j + i) 

(-IP nc20.2e(M,-^^)]^P?- V r ( \ ■ I I \ \ 

2^— p I ^ I /j+(7-i)/2(pn - J(7-i)/2(c20.2e(u,w)pr-). 



Using the formula (11.4. p. 694 in [2 ) 



m>0 



kV 



-a/2 



2z 



one finally sees that p.9p simplifies to (|c2^.20(u, w)| < 1 for almost all {u,v)) 

/ „ \ (7-l)/2 



Finally 

2^r(7) 



^(7-l)/2(P7-y^l - clg{u,v)) = i(^_l)/2(pr-^l - c2g(M,w)). 



(pr)-)- 
(2.10) 



{+)Y.{3+l)I,+-,{pr)C] {C2^.2e{u,v)w) {±zy^ p^-<~^^'\dw) = 2T 



j>0 



7 + 1 



COS^ (y-) «(7~l)/2(P''\/l -c2(,(m,u)) 



which extends to all strictly positive integer values of 7 as being zero at odd values. 
Now, let 7 be an odd positive integer, then by the virtue of (12. 6|) . it remains only 
to integrate 



(2.11) I ^ Dl[Z COshiWl - Z^)]^z=c,,Mu,v)y.f^ 

for fixed u,w e [—1, 1] and t G [0,pr]. Using Leibniz's rule, one gets: 



Dl[Z cosh{t^/l - Z^)] = ZDl[cosh{t^/l - Z^)\ + -^DY^[cosh{t^/l - Z'^)], 

10 



so that an integration by parts transforms the second integral in (|2.1ip to 

- C i^J+l[cosh(tyT^)]|^^,,^,,(„,,)„/i(^+l)/2(d«;) 



7 J -I 



Hence, (|2.1ip is ahuost similar to (|2.8p so that one gets: 

(-i)<^-^'/^/'^i?r'[cosh(tv/r^)]i^.,,,^,(„,.)„M^-+^^/^(d«') = -r (^) t^^' 

q'^'^^ j{'f+i)/2ic2ci,,2e{u, v)tq)vtidq) 



and 

(-i)^^-'^/'/'^^r'[cosh(tyT^)]|z.,,,,,,(„,.)„M^^-')/'rft" = r (^^) t^-^ 

q'^^^j{i~i)/2{c24>,2e{u,v)tq)i^t{dq). 
Also, similar computations as done before show that 

and using the derivation rule [z^^^ I^^i{z)]' — z'-'^^I^{z) ([H]) with v ~ 7/2 + j, it 
then follows that 

/■., ^7^-l. . ^ u A 07/2'^ (-lVr((7/2)+j + l) ( c%.2eiu,v) \ 
j {tlV^ Ji,+,y2{c,,Mn,v)tq)Mdq) = -2-/ - g r((7 + l)/2 + , + 1),! 2 J 

/,/2+,+iW = 2-/ d?l^r((7 + l)/2 + j)(j-l)! 2 j * ^''/^+^-^') 

(2.12) 

_ ^7/2 d ^ (-im(7/2)+j) . 

did^2-.T{{^ + l)/2 + j)j\' * V2+.W|.=ci,^,K.)/2- 

Substituting 7 + 1 by 7 — 1 , one gets 

/■^-i- r r uw^A ^2V/^-V (-ipr(7/2 + j) ft , y 

j Jh-i)/2{c2^,2e{u,v)tq)i^t{dq)= l^jj Z^ r((7 + l)/2 + i)j! (^2''^'^-^^^''' ^7/2+j-iW 
so that 

^2 



^7-l . f ( U ^ U \ 07/2-1 t« (-lVr(7/2+j) / C,^,29(".^) 
(i'Z)^ J^,.,y2ic2,Mu,v)tq)Mdq)-2^/ + + [ ^ 

(2.13) r/2+^V2+,(i)- 



Set 



J7/2+J 



(0 



E 



(-lpr((7/2)+j) /zt 



2-^/%t^r((7 + l)/2 + j-)j! 



j-2 \ J 



for 7 > and < z < 1/2, then a glance at (|2.12p and (|2.13p shows that 



(-1) 



(7-l)/2 



y" i:>j[Zcosh(a - 



(7-l)/2 



{u,v) 



7 + 1 

7 



^20,20 



7+1 



Therefore, with regard to (|2.6[) . one has for odd 7 > 1 

2^r(7) 

(pr)7 
(2.14) 



•^"1 s=l i>0 



7 + 1 



2/7^ 



V 2 y 



/ Io{prt)dt < f 



jK^ prt, 



{u,v) 



7+1 
7 



dzK^ prt, 



^20,29 



(w, u) 



which extends to all strictly positive integer values of 7 as the zero function at even 
values. Unfortunately, we did not succeed to derive any more easier expression for 
Aj{t, z) (may be this is not possible), nevertheless noting that 

V^r((7/2)+j) 



r((7 + l)/2+j) Jo ^/T^ 



1 q7/2+j-i 



dq 



for 7 > 0, then 



^7/2 q7/2-i ^(_i),r((7/2) + j) 



E' 



VT^^^r((7 + l)/2 + i)j! 



{qtzy 1^/2+ jit) dq 



1 q7/2-i 



1 



r- I n , .,7^A-/2(Vl - '2zq)dq 

Jo Vi-9Vl-2zg^^ 



1 ^7/2-1 



2^/2 V^7o 27/2 
With regard to (P?7)) . (P+P)) and (P+i)) . one gets the integrand of (fTTTU)) for all 
strictly positive integers 7 which we shall extend in the following subsection to 
strictly positive real values of 7. 

2.5. Extension to 5R(7) > 0. (|1.10p extends to all complex values of 7 lying in the 
open right half-plane in a similar way as in remark p. 194 in [Q. The extension needs 
an exponential growth of both sides of (|1.10p viewed as functions of the variable 7. 
To see this, we start with rewriting the LHS of (|1.10p as 

/' E E(j + 7)*.+7(Hc7 ic2^Mu, v)w) [^y e>(^-^)/'(rf«;), 

•^^1 s=l j>0 



then, we use the bound ([IT] p. 14) 



els>WI 



which, together the definition of i^, yields 



|v(z)| < fT^— T ^ < 



r(^. + i) ' - ^'-rH 

for z G M and positive large enough v. Thus, the LHS of (|1.10p is bounded by 

Using the bound (|2.3p . Stirling's formula and since fjS^^'^^/'^ has unit mass, then 
the LHS of (|1.10p is easily seen to be bounded by 

which is of order 0(e'^). Coming to the RHS of (fTTU)) . note that r({7+l)/2)i(^_i)/2(z) 
is boimded by e'^' (as a function of 7 for fixed 2; € M), thereby so are (|2.7|) and 
(f2l0| . Finally, since = e^^°^* e]0, 1] for aU t e [0, 1] and again T{v + l)v(z) is 
bounded by e'^', then one only needs to focus on 

r((7 + i)/2)r((7/2) + j) 
[r((7 + i)/2 + j)]2 

for all fixed j > 0, which is even bounded as a function of 7 by Stirling's formula. 
Thus, (|1.10p holds for all 7 > 0. Finally, (|1.10[) follows after integrating with respect 
to the symmetric Beta distributions , Ai'' displayed in (EJ) and the proof of (fTTOl) 
is finished. 



Remark. The function Z i~* cosh(tvl-^Z^) appeared in ( 1 ) in relation to the 
famous Bessis-Moussa-Villani's conjecture, where the author used the Dirac-Pauli 
matrices: 

then it is easy to see that 

cosh(iVl - Z-^) = tr[e*('5i+'*^'5^)]. 

Moreover, in the Lie group scope, Qi and Q2 are the generators of the B2-Weyl 
group of the compact symplectic group ([3)- 



2.6. Extension to 7 = 0. It is obvious that (|2.7p and (|2.10p have no singularity 
at 7 = 0. Hence, one only needs to deal with (|2.1ip where the only trouble comes 
from the constant term [j = 0), more precisely from r(7/2) in (note that dzA^ 
does not contain that factor since the series starts from j — 1). But, it is known 
that sin(77r/2)r(7/2) has a removable singularity at 7 = since the mirror formula 
holds 

r(z)r(i - z) = 0<z<l. 

smTTZ 

Moreover, since (sin77r/2)/7 is well defined at 7 = 0, then the integrand in (|2.14p 
may be extended to 7 = and even tends to zero as 7 does. Since sin(77r/2) 
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remains in front of the integrals after the extension, then (|2.14p tends to as 7 
does. As a result 



lim D]^{x,y)=- lim T ( 



J J [H'r-i)/2[prc24,,2e{u,v)] 



+ cos^(77r/2)i(^_i)/2[pr,/l-c2 pt'i (du)^'" (dw) 



Since 



-20 



.2e(w>«) = 



1 — u cos 20 cos 2(p — V sin 2^? sin 2(j) 



then p.7p and (|2.9p give the same contribution up to the factor cos^(77r/2). Thus, 



D^{x,y)^ hm /' r 



7 + 1 



«(7-i)/2[p'"C20,2e(M, w)]^'! {du)^i^°{dv). 



(in fact, is an even function so that ii,{\c2cj,,20{u, v)\) = ii/(c20,2e(u, «)))■ Now, it is 
an easy exercice to see that the symmetric Beta distribution /i" converges weakly as 
v — !■ —1/2 to the symmetric Bernoulli distribution (use the integral representation 
of the spherical Bessel function j^) 

ri{du) := — [5^i{du) + 5i{du)] . 

It follows from the integral representation of i[^-i)/2 that 



D^{x,y)^ hm 



cosh(prc20,20(", w)^)^^*'^ ^^^'^idz)^^^{du)ii'-"{dv) 



lim } 



1 



m=0 ^ 

fi^'^-^')/^{dz)n^'{du)fi''>{dv) 



1 + u COS 26* cos 2(j) + V sin 20 sin 20 



(prz) 



2m 



cosh(prc20,2e (w, v)z)ri{dz)ri{du)7]{dv) 

cosh[prc20,2e('«, v)]ri{du)r]{dv) 

where we used Tonelli's Theorem to exchange the order of integration coming to 
the second equality, we used the fact that the weak convergence is equivalent to 
the convergence of moments for compactly-supported distributions together with 
the boundedness of the moments of by 1 and Lebesgue's convergence Theorem 
to derive the third equality. Since 



C20,2fl(w, w) 



1 + Z2^.2e{u,v) 



^20,26 (w, u) — u cos 29 COS 20 + w sin 29 sin 20, 



and since cosh is an even function, then elementary trigonometric identities 



1 + cos(2z) 



cos^ z, 



1 — cos(2z) 



sin^ z 



yield 



[cosh{pr cos(0 + 6*)) + cosh(pr cos(0 — 6*)) + cosh(pr sin(0 + 6)) + cosh(pr sin(0 — 6))] 



which fits (|1.7p already derived in the introductory section. 
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3. Proof of Proposition (II. 3p 

Let K, m be positive integers, | • | be the Euclidean norm in and Y,n is the 
ordinary spherical harmonic of degree m written in polar coordinates as (see [TB] 
for instance) 

Y^{y) = \yr cos me, y ^ \y\e'' . 
Recall also that the generalized Bessel function is defined by (jH]) 



D^{x,y) i ^ D,ix,wy) = ^ E ^'^ 



,yw) 



iy), 



where Dk , Vk denote the Dunkl kernel and the Dunkl intertwining operator respec- 
tively ( 9 ). With the above notation, (|1.4[) reduces with p = 2 to 

D^ix,y)^lY: e(-y-^^Ioip\y\) + 2f:'-^l^YMcos{m- 

weB2 j=i 

Consequently, by linearity of Vfc (PJ), one has 



(3.1) Df'{x,y) = V, [/o(p| • I)] iy) + 2j2Vk 



hM ■ I) 



Y 



4j 



(y) cos(4j(/)). 



Clearly we can expand p'"^''" (cos (/)) as a sum of cosmcj), < m < j so that, 



p'j'''°{cosA(j)) — E bjnj cos{Am(j)) 



m— 



(the coefficients bmj may not be simple in general, although it has nice formula 
in the case of li — Iq, that is, for Gegenbauer polynomials). Substituting this 
expression into taken with p = 2 and changing the summation's order, we get 



oc / oo 



D^{x,y) = C2.fc f -j E bm.jhj+y{pr)p'l''°{cos{'ie)) | cos(4m0). 

^ ^ ' m—O \.j—rn 



Comparing the last equation with (j3.ip . we obtain upon using the orthogonality of 
cos Arriip that 



(3.2) 



Vk 



hrn{p\ ■ I) 



Ya, 



p {p\y\)' ^ 



Using the expansion of the modified Bessel function \1.2\ to expand Iij+-y in 
the RHS of becomes 



E E 



1 



C2,fe 



oo oc 
9=0 i=0 



- g!r(4j + 7 + 9+1) 



2 



2q+4j 



p'/^'"(cos(40)) 



E 



E 

q+2j=K 



'g!r(4j + 4m + 7 + g+l) V 2 



g!r(4j + 4m + + 7 + 1) V 2 



2g+4j+4m 



p;-'^(cos(40)) 



2K+4m 



p'liticosm) 



15 



Similar computations transform the Bessel function in the LHS of p.2p to 

hra{p\-\) _ 1 1 

(p| • 1)4™ 24™ ^ K!r(4m + K + 1) V 2 y ^ 
Comparing the coefficients of the obtained expressions, one finally gets 

vu [I ■ rn^] {y) ^ ,^^-^ucodm 

K!r(4m + K + f ) ^ (K-2j)!r(4m + 2j + «; + 7+l)^^+™^^ " 

^ 2 J ^ K, 

and Proposition (|f .3p is proved. 
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